Universal and nonuniversal allometric scaling behaviors in the visibility graphs of 

world stock market indices 
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The investigations of financial markets from a complex network perspective have unveiled many 
phenomenological properties, in which the majority of these studies map the financial markets into 
one complex network. In this work, we investigate 30 world stock market indices through their 
visibility graphs by adopting the visibility algorithm to convert each single stock index into one 
visibility graph. A universal allometric scaling law is uncovered in the minimal spanning trees, 
whose scaling exponent is independent of the stock market and the length of the stock index. 
In contrast, the maximal spanning trees and the random spanning trees do not exhibit universal 
allometric scaling behaviors. There are marked discrepancies in the allometric scaling behaviors 
between the stock indices and the Brownian motions. Using surrogate time series, we find that 
these discrepancies are caused by the fat-tailedness of the return distribution, the nonlinear long- 
term correlation, and a coupling effect between these two influence factors. 

PACS numbers: 05.45.Tp, 05.40.-a, 05.45.Df, 89.75.Da, 89.65.Gh 



I. INTRODUCTION 

Econophysics is an interdisciplinary field which adopts 
ideas, tools and theories from statistical mechanics, non- 
linear science, complexity science and applied math- 
ematics to understand the emerging complexity and 
self-organized macroscopic behaviors of economic sys- 
tems, with special interest paid to financial markets 
EL BIS BL Hi ■ Econophysists are particularly interested in 
unveiling different universal behaviors in financial mar- 
kets [H, 0, M . [si [TpL ITlT ] , following the phenomenological 
framework [12 . 1 1 31 ] . 

In recent years, com plex network theory has witnessed 
a flourishing progress [3 G3, QJ, G3, GJ, 0- It is nat- 
ural that a wealth of studies have been carried out from 
a complex network perspective. The current economic 
crisis calls for a deeper understanding of the dynamics of 
economic activities on the global economic network [20l ]. 
The studies in this field can be classified into two types 
based on how the network is constructed. The studies of 
the first type deal with many time series to form a com- 
plex network with each node standing for a time series 
and the weight of a link between two nodes character- 
ized by the correlation coefficient of the two time series 
2l| . [22. |23| or by the distance between the two time series 

Concerning the studies of the second type, different 



mapping methods have been proposed to convert time 
series into different kinds of networks, including cycle 
networks based on the local extrema and their distance 
in t he p hase space [27], [28| , segment correlation networks 
[H,[l0|, nearest neighbor network s l3ll . n -tuple networks 
based on the fluctuation patterns [32|, [33| , the visibility of 
nodes [34[ , space state networks based on conformational 
fluctuations [35l. b in transition networks [36| , and recur- 
rence networks [37], [H| . The visibility algorithm has been 
diversely used to investigate stock market indices [39| . 
human strive intervals [40 ], occurrence of hurricanes in 
the United States [4l[, foreign exchange rates j^, and 
energy dissipation rates in three-dimensional fully devel- 
oped turbulence [43j . 



'Electronic address: wxzhou@ccust.cdu.cn 



In this work, we study the allometric scaling behavior 
of spanning trees extracted from the visibility graphs of 
30 stock market indices all over the world. Both uni- 
versal and nonuniversal scaling behaviors are reported, 
which is reminiscent of the universal scaling behavior of 
the weighted world trade networks [44|. The paper is 
organized as follows. In Sec. UH we describe briefly the 
methodology adopted. Section [TTT1 presents the empirical 
findings. We study the impact of the length of the stock 
indices on the allometric behaviors in Sec. IIVI Section 
IVl explores the driving factors that cause the different 
behaviors between the stock indices and the Brownian 
motions. And Section IVll concludes. 
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II. METHODOLOGY 



C. Allometric scaling 



Here we briefly explain the methodology used in this 
study. For each stock market index, a unique visibility 
graph is constructed. A maximal spanning tree (MaxST), 
a minimal spanning tree (MinST), and 100 random span- 
ning trees (RanSTs) are extracted from the constructed 
visibility graph. For each spanning tree, an allometric 
scaling analysis is carried out and a scaling exponent is 
determined. 



A. Construction of visibility graph 

Consider the price series p(t) of a stock market index 
with a length of N. We can transform the time series 
into c omp lex networks by applying the visibility algo- 
rithm [34[ . Each data point in the time series is regarded 
as a node in the complex network, and an edge is drew 
connecting two nodes according to the rule that the two 
corresponding data points can see each other in the di- 
agram of the time series. Mathematically, two arbitrary 
data points p(ti) and p(tj) have visibility if any other 
data point p(tk) located between them fulfills 



P(tj) - P(tk) ^ p(tj) ~p{tj) 



U - t* 



> 



1 7 ti 



(1) 



We assign the average acceleration of price movement, 
defined as [hxp(tj) — lnp(i,)]/(ij — ti), as the weight of 
the edge [ti,tj], where U < tj. 



Allometric scaling laws are ubiquitous in complex 
systems evolving on complex networks, such as the 
metabolism of organisms and ecosystems river networks 
[H, H| 53, JM 51 HI, the food webs [H, the world 
trade webs |44j . the world investment networks [52| , and 
so forth. The original model of the allometric scaling on 
a spanning tree was developed by Banavar, Maritan, and 
Rinaldo [49j | . The node with the maximum degree is con- 
sidered as the root of a spanning tree. Each node of a 
spanning tree is assigned a number 1, and two values Ai 
and Ci are defined for each node i in an iterative manner 
as follows: 



+ 1 and d = Y^Cj + Ai 



(2) 



in which j stands for all the nodes linked from i [49|| . 
The allometric scaling relation is then highlighted by the 
power law relation between Ci and Af. 



C~A V , 



(3) 



where the leaf nodes with A = C = 1 should be excluded 
from the estimation of the scaling exponent r\ (5l| . 

Any spanning tree can range in principle between two 
extremes, that is, the chain-like trees and the star-like 
trees. For chain-like trees we have r\ = 2~, while for 
start-like trees we have ?; = 1 + . Therefore, 1 < rj < 2 for 
all spanning trees. It should be note that not all trees 
exhibit such an allometric scaling behavior, for instance 
the classic Cayley trees [53[ . 



B. Construction of spanning trees 

We extract three different kinds of spanning trees 
(MaxST, MinST, RanST) from each visibility graph. The 
algorithms for building MaxST, MinST, and RanST are 
given in the following. 

MaxST: First, select an edge with a maximum weight 
as the first edge of the MaxST. Second, select a 
new edge with a maximum weight among the edges 
which connect the tree and make sure that no loop 
is introduced. Third, repeat the second step till all 
nodes are added into the tree. 

MinST: First, select an edge with a minimum weight 
as the first edge of the MinST. Second, select a 
new edge with a minimum weight among the edges 
which connect the tree and make sure that no loop 
is introduced. Third, repeat the second step till all 
nodes are added into the tree. 

RanST: First, arbitrarily select an edge as the first edge 
of the RanST. Second, randomly select a new edge 
among the edges which connect the tree but make 
sure that no loop is introduced. Third, repeat the 
second step till all nodes are added into the tree. 



III. UNIVERSAL SCALING IN THE 
VISIBILITY GRAPHS CONSTRUCTED FROM 
STOCK MARKET INDICES 

The data sets we analyzed contain 30 stock market 
indices all over the world, which are retrieved from Ya- 
hoo! Finance at http://financc.yahoo.com A list of 



the index names is given below together with the ab- 
breviations, countries or areas and starting dates of the 
time series used for analysis in the ensuing parenthe- 
ses: Amsterdam Exchange Index (AEX, Netherlands, 3 
January 2000), ATX Vienna (ATX, Austria, 11 Novem- 
ber 1992), Euronext BEL-20 (BFX, Belgium, 11 Febru- 
ary 2005), BSE Sensex (BSESN, India, 1 July 1997), 
Ibovespa (BVSP, Brazil, 27 April 1993), CMA GENL In- 
dex (CMA, Egypt, 26 May 2003), Dow Jones Industrial 
Average (DJIA, USA, 1 October 1928), CAC 40 Index 
(FCHI, France, 1 March 1990), FTSE 100 Index (FTSE, 
UK, 2 April 1984), DAX Index (GDAXI, Germany, 26 
November 1990), S&P/TSX Composite Index (GSPTSE, 
Canada, 3 January 2000), Hang Seng Index (HSI, Hong 
Kong, 31 December 1996), Jakarta Composite Index 
(JKSE, Indonesia, 1 July 1997), FTSE Bursa Malasia 
KLCI Index (KLSE, Malaysia, 3 December 1993), Ko- 
rea Composite Stock Price Index (KOSPI, Korea, 1 July 
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1997), Mcrval Buenos Aires (MERV, Argentina, 8 Octo- 
ber 1996), MIBTEL Index (MIBTEL, Italy, 3 January 

2000) , MXX IPC (MXX, Mexico, 1 November 1991), 
NIKKEI 225 (N225, Japan, 4 January 1984), NASDAQ 
Composite (NASDAQ, USA, 2 February 1971), NZX 
50 Index (Gross) (NZ50, New Zealand, 30 April 2004), 
OMXS All Share Index (OMXSPI, Sweden, 8 January 

2001) , Oslo Exchange All Share Index (OSEAX, Nor- 
way, 7 February 2001), IGBM (SMSI, Spain, 2 January 

2002) , Standard and Poor's 500 Index (S&P500, USA, 
3 January 1950), Shanghai Stock Exchange Composite 
Index (SSEC, China, 4 January 2000), Swiss Market In- 
dex (SMI, Switzerland, 9 November 1990), Straits Times 
Index (STI, Singapore, 28 December 1987), Tel Aviv TA- 
100 Index (TA100, Israel, 1 July 1997), and Taiwan Stock 
Exchange Corporation Weighted Index (TWII, Taiwan, 
2 July 1997). The ending dates of all the indices are 25 
August 2009. 

For each stock index, a visibility graph is constructed 
and its maximum spanning tree (MaxST) and minimum 
spanning tree (MinST) are determined uniquely. In ad- 
dition, 100 random spanning trees (RanSTs) are also de- 
rived from the visibility graph. For each tree, an allomct- 
ric analysis is carried out and the two sequences of A and 
C are calculated. Figure [1] shows the allometric scaling 
behaviors of the MaxST, the MinST and a randomly se- 
lected RanST associated with the FSTE 100 Index. Nice 
power-law relationships are observed between A and C. 
A linear regression to the data finds that ?/ = 1.271±0.002 
for the MaxST, rj = 1.264 ± 0.002 for the MinST, and 
r/ = 1.308 ± 0.002 for the RanST, respectively. 




A 

FIG. 1: Allometric scaling behavior of spanning trees ex- 
tracted from the visibility graph of the FTSE 100 Index. The 
data points for MinST and RanST are transformed vertically 
by a factor of 10 and 0.1 for better visibility. The solid lines 
are the best power-law fits. 

We find that all the spanning trees exhibit excellent 
allometric scaling behaviors for all the indices. The cor- 
responding power-law exponents for different indices are 
reported in table [J We find that the allometric scaling 
exponents of the MaxSTs for different indices are very 
close to each other, centered around 1.281 ± 0.009 (mean 



± std). The situation for the MinSTs is the same as 
the MaxSTs, in which the power-law exponents fluctuate 
slightly around 1.274 ±0.011. In contrast, the power-law 
exponents of the RanSTs are basically larger than 1.3 
except for DJIA, NASDAQ and S&P 500, which are sig- 
nificantly larger than the exponents of both the MaxSTs 
and the MinSTs. We infer that the appearance of smaller 
values of r? RanS T for DJIA, NASDAQ, and S&P 500 is due 
to the fact that these three indices have much more data 
points compared with other indices. 



TABLE I: Allometric scaling exponents rj of spanning trees 
(MaxST, MinST, and RanST) for different indices. 



stock index 


??MaxST 


^MinST 


TjRanST 


AEX 


1.276(3) 


1.278(3) 


1.329(3) 


ATX 


1.282(3) 


1.270(3) 


1.312(2) 


BFX 


1.270(5) 


1.269(4) 


1.349(4) 


BSESN 


1.288(3) 


1.281(3) 


1.310(3) 


BVSP 


1.299(3) 


1.290(3) 


1.315(2) 


CMA 


1.277(12) 1.250(8) 


1.315(7) 


DJIA 


1.284(1) 


1.267(1) 


1.274(1) 


FCHI 


1.280(2) 


1.269(2) 


1.313(2) 


FTSE 


1.271(2) 


1.264(2) 


1.308(2) 


GDAXI 


1.284(2) 


1.270(2) 


1.312(2) 


GSPTSE 


1.273(3) 


1.260(4) 


1.333(3) 


HSI 


1.276(2) 


1.276(2) 


1.304(2) 


JKSE 


1.284(3) 


1.280(3) 


1.309(3) 


KLSE 


1.276(3) 


1.284(3) 


1.318(2) 


KOSPI 


1.279(3) 


1.292(3) 


1.307(3) 


MERV 


1.278(3) 


1.282(3) 


1.318(2) 


MIBTEL 


1.274(3) 


1.274(3) 


1.327(3) 


MXX 


1.281(3) 


1.280(3) 


1.312(2) 


N225 


1.276(2) 


1.271(2) 


1.308(2) 


NASDAQ 


1.308(2) 


1.267(2) 


1.271(2) 


NZ50 


1.260(4) 


1.276(4) 


1.348(4) 


OMXSPI 


1.278(3) 


1.283(4) 


1.313(3) 


OSEAX 


1.278(3) 


1.266(4) 


1.309(3) 


SMSI 


1.297(5) 


1.294(5) 


1.345(5) 


SP500 


1.280(1) 


1.256(1) 


1.282(1) 


SSEC 


1.285(4) 


1.264(4) 


1.313(3) 


SMI 


1.279(2) 


1.267(2) 


1.313(2) 


STI 


1.273(2) 


1.276(2) 


1.317(2) 


TA100 


1.290(3) 


1.286(3) 


1.324(3) 


TWII 


1.291(3) 


1.285(3) 


1.325(3) 



IV. FINITE-SIZE EFFECT 

Table Q] shows that the allometric scaling exponents 
for a given type of spanning tree are very close to each 
other. However, there are still fluctuations around the 
corresponding average values. It is possible that the ex- 
ponent is dependent on the length of the index time se- 
ries. It is thus necessary to further investigate if there 
is a finite-size effect, which is crucial to the validation of 
universality. 



4 



A. The case of MaxST and MinST 

Figure [5] illustrates the dependence of r^MaxST and 
7?MinST as a function of the stock index length L. No 
evident trend is identified by eye-balling. We fit the data 
for MaxSTs and MinSTs using a linear model 

?7MaxST\MinST = a + bL. (4) 

For the MaxSTs, a = 1.279 and b = 3.949 x 1CT 7 and 
the corresponding p- values from the Student's t-test are 
and 0.355, respectively. For the MinSTs, a = 1.278 
and b = —8.315 x 10 and the corresponding p-values 
are and 0.082, respectively. If we regress the exponent 
against InL, the p- value of b is 0.322 for the MaxSTs and 
0.314 for the MinSTs. It is clear that the coefficient b is 
identical to and the exponents ^MaxST and ^MinST are 
independent of L. 



S=~ 1.28 
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FIG. 2: (Color online) Dependence of the allometric scaling 
exponents ryiviaxST and ryuiinST on the length L of stock indices. 
No evident finite-size effect is observed. 



B. The case of RanST 

For the RanSTs extracted from the visibility graphs 
of stock indices, the three allometric scaling exponents 
?7Ra„ST of the US market indices (DJIA, S&P 500, NAS- 
DAQ) are significantly less than the others, which is a 
signal of the possible presence of finite-size effect. Figure 
[3] plots ?7RanST as a function of L for all the indices. It 
is evident that ?7RanST decreases logarithmically with the 
increase of index length: 

?7RanST = a + bh\L, (5) 

where a = 1.471 and b = —0.019, both of which are 
significantly different from with the p-values less than 
IO" 6 . 



1.38 
1.36 
H 134 

CO 

§1.32 
K- 1.3 
1.28 



2 3 4 5 

10 10 10 10 

L 

FIG. 3: Dependence of the allometric scaling exponent r;R an sT 
on the length L of stock indices. An evident finite-size effect 
is observed. 

C. Numerical tests 

In order to further investigate the impact of the time 
series length on the allometric scaling exponents of the 
three types of spanning trees, we design and conduct two 
numerical tests. We take subseries with different lengths 
from the longest index DJIA to perform the allometric 
scaling analysis. For each length L, 100 subseries are 
randomly extracted and the corresponding scaling expo- 
nents are calculated from their MaxSTs, MinSTs, and 
RanSTs. Figure 0] illustrates the dependence of the av- 
eraged exponent rj with respect to the length L for the 
DJIA subseries. A linear regression of r^MinST against L 
using Eq. g]) gives that a = 1.261 and b = 4.093 x 10 -8 
with the p- values being and 0.16, respectively. It means 
that ?7MinST is independent of L. For the MaxSTs, we find 
that ?7MaxST increases with L linearly 

WaxST = 1-276 + 4.555 x 10- 7 L, (6) 

where the linear coefficients are statistically significantly 
different from zero with both the p-values being zero. In 
contrast, ryR an ST exhibits an evident decreasing trend. A 
linear regression gives a — 1.327 and b = —2.581 x 10~ 6 , 
whose p-valucs are both nulls. These observations are 
consistent with the results in Tabled] In addition, we find 
that TyMinST < ?7MaxST for all L's, while ?7 RanST becomes 
less than ?7MaxST when L > 16500 and less than ?7MinST 
for much greater L. 

For comparison, we synthesize Brownian motions with 
different lengths. For each length, 100 time series are 
generated and the allometric scaling exponents are de- 
termined. The results are also illustrated in Fig. [4] The 
scaling exponents ryMaxST and r^MinST are found to be in- 
dependent of L since linear regressions show that the p- 
values of the corresponding b's are 0.17 and 0.06, respec- 
tively. On the other hand, the scaling exponent ?7MinST 
exhibits a decreasing trend for small L and then reaches 
a constant. For all L's, 7/MinST ~ ?7MaxST < ??RanST for 
the Brownian motions. 
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FIG. 4: (Color online) Dependence of the averaged exponent 
77 with respect to the length L for the DJIA and the Brownian 
motions. 



We stress that there are marked discrepancies between 
the scaling exponents of the MinSTs (and the MaxSTs 
as well) of the DJIA and Brownian motions. It calls 
for an investigation of influence factors on the allometric 
behaviors of the MinSTs and MaxSTs of the stock market 
indices. Note that a comparison between the behaviors 
of the stock indices and the Brownian motions is not 
out of blue. Indeed, stock prices are assumed to follow 
a geometric Brownian motion in the celebrated Black- 
Scholes model of option pricing [54| . 



V. INFLUENCE FACTORS ON THE 
ALLOMETRIC SCALING BEHAVIOR 

To understand the marked discrepancies between the 
allometric behaviors of the stock market indices and the 
Brownian motions, further numerical experiments are 
needed. Comparing with Brownian motions, we find 
there are three potential factors that may have influ- 
ence on the behavior of any time series, that is, the fat- 
tailcdness of the probability distribution, the linear long- 
term correlation, and the nonlinear long-term correlation 
[55l ]. It is well established that there is no linear long- 
term correlation in the returns of stock indices. Hence, 
two factors remain for further investigations. For com- 
pleteness, we first study the impact of linear long-term 
correlations based on fractional Brownian motions. 



square linear fit to the data for each spanning tree yields 
an estimate of the power-law exponent, which results in 
r) = 1.233 ± 0.004 for the MaxST, 77 = 1.233 ± 0.004 for 
the MinST, and 77 = 1.314 ± 0.004 for the RanST. Note 
that the node (1,1) is excluded in the implementation of 
fitting. The plots of A with respect to C for the other 
Hurst indexes share the similar pattern as the plots for 
H = 0.5, but the power-law exponents arc different. 

Figure [5] shows the allometric scaling exponent 77 as a 
function of the Hurst index H for three different spanning 
trees. For the MaxST and the MinST, the two curves 
almost overlap onto a same curve. Furthermore, both 
curves show a linear decreasing trend 



rj = a + bH, 



(7) 



where a = 1.238 and b = -0.008 for the MaxSTs and 
a = 1.238 and b = -0.007 for the MinSTs with all the p- 
valucs less than 0.01. The variation of these two scaling 
exponents is actually very slight and all the exponents 
are embedded in the interval [1.23, 1.24]. In contrast, the 
power-law scaling exponent 77 for the RanSTs increases 
with the Hurst index H, which can also be modeled by 
Eq.([7]). A linear regression finds that a = 1.301 and 
b = 0.027 with the associated p-values being zero. 
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FIG. 5: Plots of the dependence of allometric scaling exponent 
77 on the Hurst index H for three spanning trees. 



B. Origin of the difference of allometric scaling 
exponents between financial series and Brownian 
motions 



A. Fractional Brownian motions 

We generate fractional Brownian motions through 
wavelet transform with the Hurst indexes varying from 
0.05 to 0.95. We obtain 100 realizations for each Hurst 
index, and each realization has 5000 data points. 

In the case of Brwonian motions where the Hurst index 
H = 0.5, remarkable power-law behaviors are observed 
between A and C for the three spanning trees. The least 



Compared with the Brownian motions, the financial 
series exhibit fat-tailed PDF and nonlinearity as well [55[ . 
We design several tests to verify the contribution of the 
two factors on the allometric scaling behavior. In the 
tests, different surrogate data are used. There is no linear 
long-term correlation in the surrogate data. We use the 
FTSE index as an example, and all the surrogate time 
series have the same length as the FTSE data. 

The surrogates of the first type (termed "Surr 1") have 
the same probability distribution of returns as the FTSE 
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index, but without any nonlincarity. We can shuffle the 
original the corresponding return series r(t) to remove the 
nonlincarity, and conduct cumulative summation to re- 
construct the prices. Alternatively, we can determine the 
empirical distribution of the FTSE returns and generate 
surrogate return series (56[. We find that the surrogate 
time series obtained from either methods give the same 
scaling exponents. 

The surrogates of the second type (termed "Surr 2") 
preserve the nonlinearity of the FTSE data but with 
Gaussian return distributions. We generate a Brownian 
motion, rearrange its increment series to ensure that the 
resultin g se ries has the same rank ordering as the FTSE 
returns |55l.[57l.[58l]. 

The surrogates of the third type (termed "Surr 3") 
preserve both the nonlinearity and the probability dis- 
tribution of the original FTSE returns. We synthesize a 
surrogate return series based on the empirical distribu- 
tion of the FTSE returns and then rearrange the rank 
orders of the data to introduce nonlinear correlations. 

For each test, 100 surrogates are generated, and the 
allomctric scaling exponents are determined for the cor- 
responding MaxSTs, MinSTs, and RanSTs. The results 
are presented in Table [Til which are compared with the 
FTSE index and the Brownian motions. First of all, the 
scaling exponents of the random spanning trees are close 
to each other for the five types of data. For the MaxSTs 
and MinSTs, we have 

VBm < '/Suit 1, ?/Surr 2 < ?7Surr 3 = ??FTSE- (8) 

Therefore, both the nonlinearity and the fat-tailedness 
have influence on the allometric scaling. However, only 
one factor is not capable of explaining the difference of rj 
between the FTSE index and the Brownian motions. A 
coupling effect between the two factors also has contri- 
bution. 



VI. CONCLUSION 



In summary, we have investigated the allometric scal- 
ing behaviors of the maximal spanning trees, the minimal 
spanning trees and the random spanning trees of the visi- 
bility graphs constructed from 30 worldwide stock indices 
with different sizes. All the spanning trees exhibit nice 
allometric scaling behaviors. We found that the aver- 
age scaling exponent is ?7MaxST = 1-281 ± 0.009 for the 
maximal spanning trees and r/MinST = 1-274 ± 0.011 for 
the minimal spanning trees. Using subscries extracted 
from the longest index series of the DJIA, we found that 
the exponent for the minimal spanning trees is indepen- 
dent of the length of the index, while the exponent for the 
maximal spanning trees exhibits a weak increasing trend. 
In contrast, the exponent for the random spanning trees 
decreases fast with the increase of the length. Therefore, 
the visibility graphs of the world stock market indices ex- 
hibit both universal and nonuniversal allomctric scaling 
behaviors. 

Numerical simulations showed that Brownian motions 
with different sizes exhibit different allometric behaviors 
in the associated maximal, minimal, and random span- 
ning trees. Surrogate time series with different features 
were synthesized to study the origin of the discrepancy 
between Brownian motions and stock indices. We con- 
clude that the fat-tailedness of the return probability 
distribution, the nonlinear long-term correlation, and a 
coupling effect between them can explain exactly the dis- 
crepancy. 



TABLE II: Determining the influence of the fat-tailedness in 
the probability and the nonlinear long-term correlation on the 
allometric scaling behavior of the FTSE index. 



Test MaxST MinST RanST 

FTSE 1.271 ±0.002 1.264 ±0.002 1.308 ±0.002 

Surr 1 1.263 ±0.004 1.250 ± 0.005 1.306 ±0.005 

Surr 2 1.255 ±0.002 1.255 ± 0.002 1.315 ±0.004 

Surr 3 1.272 ±0.002 1.262 ± 0.003 1.310 ±0.004 

Bm 1.241 ±0.002 1.243 ± 0.002 1.303 ±0.004 
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